DIFFUSIVE WAVELETS ON GROUPS AND HOMOGENEOUS SPACES 



SVEND EBERT AND JENS WIRTH 



Abstract. The aim of this exposition is to explain basic ideas behind the concept of diffusive 
wavelets on spheres in the language of representation theory of Lie groups and within the 
framework of the group Fourier transform given by Peter- Weyl decomposition of L^{5^) for a 
compact Lie group 'S . 

After developing a general concept for compact groups and their homogeneous spaces we 
give concrete examples for tori — which reflect the situation on R" — and for spheres and S''. 
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1. Introduction 

Within the past decades wavelets and associated wavelet transforms have been intensively used in 
both applied and pure mathematics. They and the related multi-scale analysis provide essential 
tools to describe, analyse and modify signals, images or in rather abstract concepts functions, 
function spaces and associated operators. We refer to the seminal expositions of Meyer [21] and 
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Coifman-Meyer (2- The most influential constructions in M" are due to Haar [T^, Grossmann- 
Morlet HI] and Daubechies 0. 

The constructions of wavelet transforms can be entirely based on an abstract group theoretic and 
representation theoretic approach. An abstract exposition of this can be found in Kisil [20 , for 
the particular situation of the Lorentz group SO(n + 1, 1) acting on the sphere S" an associated 
associated wavelet constructions was carried out by Antoine and Vandergheynst [I] , [2] ; see also 

m, m- 

We will shortly recall the general group theoretic approach. Let be a (pseudo-) Riemannian 
manifold. A wavelet transform in L^[Ai) is defined in terms of an irreducible unitary represen- 
tation U of Lie group 

U :^ ^ C{L'^{M)). (LI) 
A non-zero vector ^ 6 L'^{M) is an admissible wavelet if 

/ \{f,U{g)^)L2^M)\'dg <oo (L2) 

for all / e L'^{A4). The integral is taken with respect to (any) Haar measure on Equation 
(|1.2I) implies that the associated wavelet transform 

^/(5)-(/,t/(5)*)L^(A<) (1-3) 

is a bounded operator W : L^{A4) L^(J^), while the irreducibility of U ensures invertibility of 
W on its range. To see this, we take an / G ker^'. Then / ± span{?7(f/)4' : g G Since U 
is unitary, the orthocomplement of an invariant subspace is invariant and by irreducibility of U 
this implies / = 0. 

This approach has several disadvantages, which becomes apparent in particular for spheres. 
Because L^(A^) is infinite dimensional, no compact group admits an irreducible unitary repre- 
sentation of this form. However, compact groups seem natural at least in the situation where 
M itself is a homogeneous space of a compact group. A related approach was incorporated for 
spheres by Antoine and Vandergheynst. They used quasi-regular representations of the Lorentz 
group SO(n -f- 1,1) on the sphere S" ~ SO(n + 1, 1)/S0(n, 1). To overcome the the lack of 
square-integrability they restricted considerations to so-called admissible sections of the Lorentz 
group SO(n + 1, 1). See [TU], [TT] for a detailed exposition on that. 

An alternative approach to wavelet transforms on spheres was followed by Coifman-Maggioni 
in [B] or by Freeden [T2]. Classical wavelet theory on M" is based on the group generated by 
translations and dilations. It is evident what translations are rotations on a sphere (seen as 
homogeneous space of the rotation group), but there is no canonical choice for dilations. The 
key idea of diffusive wavelets is to generate dilations from a diffusive semigroup, e.g., from time- 
evolution of solutions to a heat equation on the homogeneous space. The use of compact groups 
is advantageous in view of the availability of powerful tools like Peter- Weyl theorem and available 
classifications of irreducible representations. A related concept based on spectral calculus of the 
Laplacian on closed manifolds was proposed by Geller, I15j . 

Discrete wavelet transforms in such a setting were discussed in [B] and [S] , where heat evolution is 
combined with a orthogonalisation procedure to model a multi-resolution analysis within L^(§'^). 
Our approach to continuous wavelet transforms given in the present paper follows the discussion 
of |12j for spheres, but replacing the heat flow by a more general approximate convolution identity. 

In the following sections we present a representation theory based approach to understand dif- 
fusive wavelet transforms on compact Lie groups and their homogeneous spaces. The wavelet 
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transforms itself are constructed in terms of kernels of an approximate convolution identity in 
L^C^), ^ a compact Lie group, or in L^(^^/^) with Jif <'^^ a closed subgroup of 

The main structure of the paper is as follows. In Section[2]we recall some of the main ingredients 
from representation theory, in particular Peter- Weyl theorem and the resulting non-commutative 
Fourier transform on compact groups. For a more detailed introductory account on the analytical 
aspects and relations to operator theory and partial differential equations of the group Fourier 
transform we refer to [22] ■ 

Section [3] sketches the main ideas of the approach and discusses the construction of diffusive 
wavelet transforms on arbitrary compact groups. Starting point is the notion of diffusive ap- 
proximate identities given in Definition 13. 1[ for which the heat kernel provides one interesting 
and simple example. Later all constructions will be explicitly discussed for this particular case. 
The concept of diffusive approximate identity is closely related to the scaling function within the 
classical discrete wavelet approach. 

The constructed wavelet transform of / G L^(^^) has the usual structure of an inner product 
between / and a translated wavelet i/^p, p being a scaling parameter, and also the corresponding 
inversion formula has the familiar form of a superposition of translations of ^pp integrated over 
all scales, see formulas p. 21) and p.3|) . 

A projection methods transfers these wavelet transforms to homogeneous spaces of by identi- 
fying them as quotients j M' for a suitable subgroup M' < ^. On the other hand, the wavelet 
transform can be intrinsically defined on the homogeneous space itself. Section 0] provides two 
approaches to wavelet transforms on homogeneous spaces, one giving functions parametrised by 
points in the space and one parametrised by elements of the group. They extend existing con- 
structions on spheres due to [12] or |3] . As one essential ingredient for our approach we include a 
discussion of several convolution-like products on homogeneous spaces and their relation to the 
group Fourier transform. 

Finally, Section [S] contains concrete examples of transforms arising from our methods. In partic- 
ular, we discuss the unit circle, the group of unit quaternions and spheres as homogeneous spaces 
of SO(n +1). Furthermore, real projective space can be considered analogously to spheres. 



2. Prerequisites: Representation theory of compact Lie groups 

We start by recollecting some basic notation. Let be a compact Lie group. A unitary represen- 
tation of is a continuous group homomorphism tt : ^ ^ U(c?7r) of into the group of unitary 
matrices of a certain dimension dj^. Such a representation is irreducible if TT(g)A — Air^g) for all 
g G 5^ and some A G C'^^^'^^ implies A = cl is a multiple of the identity. Equivalently, C*" does 
not have non-trivial 7r-invariant subspaces V C C*" with TT{g)V C V for all 5 G We call two 
representations tti and tt2 equivalent, if there exists an invertible matrix A with Tri{g)A ~ ATT2{g) 
for ah g e'^. 

Let denote the set of all equivalence classes of irreducible representations. Then this set 
parameterises an orthogonal decomposition of L'^('i^), where d/ig^ is the normalised Haar measure 
on ^. 

Theorem 2.1 (Peter- Wevl.|24[ Chapter 1.4.3]). Lef^^ be a compact Lie group. Then the follow- 
ing statements are true. 
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(1) Denote = {g ^ iv&ce{'K{g)A) : A € C'-x'*-}. Then the Hilhert space L'^{'^) decom- 
poses into the orthogonal direct sum 

L'm^^H^. (2.1) 

(2) For each irreducible representation it the orthogonal projection L^{^) — >■ is given 
by 

^>-^d^ (f){h)xn{h~'^g)diJ,c^{h) = d^(f)*Xn (2.2) 
in terms of the character x-^{9) = trace 7r(g') of the representation. 

We will denote the matrix A in cp * Xtt = tia.ce{Tr{g)A) as Fourier coefficient 0(7r) of at the 
irreducible representation tt. A short calculation yields 

k^)= I (P{g)TT*{g)dti^{9) (2.3) 

and therefore 

<i){g) = J2d^ trace(7r(5)(^(7r)). (2.4) 

If we denote by ||>l||//5 = trace(>l*A) the Frobenius or Hilbert- Schmidt norm of a matrix A, 
then the following version of Parseval identity holds true. 

Corollary 2.2 (Parseval- identity) . Let <p G L'^{^). Then the matrix-valued Fourier coefficients 
^(tt) e C''-^''- satisfy 

\\<i>f = Y.dM'^)\\ls. (2.5) 



The group structure of ^ defines an involution of functions on We denote it as 4>{9) = 4'{9~^)- 
On the level of Fourier coefficients this corresponds to taking adjoints of the matrix- valued Fourier 
coefficients. 

On the group ^ one defines the convolution of two integrable functions (f),'ip G L^{^) as 

<t>*^p{9)= [ cj>{h)i^{h-^g)Aii^{h). (2.6) 

Since (f)*il) G L^{^), the Fourier coefficients are well-defined and they satisfy 
Lemma 2.3 (Convolution theorem on '^). Let (p,tp G i^(?f) then (j)*ip G L^($f) and 

f^{Tr) = ■0(7r)(^(7r). (2.7) 

Proof. We consider only the convolution formula. As short calculation yields 

0*i/;(7r) = jj (j){h)'4){h~^g)-K*{g)diJ.<^{h)dii<^{g) 

ct>{h)i,{h-^g)n*{h-^g)Tr*{h)dn^{h)dn^{g) = i>{Tr)${w). 



-If 



□ 
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The group structure gives rise to left and right translations, Tg : <^(ffO and : (j) ^ (t){-g ^) 
of functions on the group. For completeness, we remind the reader to the formulae 

= <^(7r)7r(g), and f^(7r) = ^*(g)^(7r). (2.8) 

They are a direct consequence of the definition of the group Fourier transform (|2.3p . (|2.4p . The 
assignments g H> Tg and g ^ are usually referred to as the left respectively right regular 
representations of in L^{'^). Both representations are unitary, T* = Tg-i and (T^)* — T^^ ^\ 

We need one more ingredient for later considerations. The Laplace-Beltrami operator Ag? on the 
group ?f is bi- invariant, i.e. it commutes with all Tg and . Therefore, all of its eigenspaces are 
bi- invariant subspaces of L^(^^). As T-Lt^ are minimal bi-invariant subspaces, each of them has to 
be eigenspace of Ac^ and we denote the corresponding eigenvalue as — A^. Hence, we obtain 

l^cgcj) = - ^ d^Xl trace(7r(.g)(/)(7r)) (2.9) 

and the solution to the heat equation 

(9t-A^)u = 0, u{{),-)^4> (2.10) 
is given as convolution with the heat kernel pt{g) as u{t, •) — (j) *pt with 

Pt(^) = e-*^'l, pt{g) = d.e-*^'x^(ff)- (2.11) 

7re# 

We remark that in particular (j) * pt ^ (j) ior all G L^'(^^), p G [1, oo). This fact will be crucial 
for the proof of inversion formula later on. 

3. Diffusive wavelets on a compact Lie group 

3.1. General philosophy. Let pt G L^{'^) be an approximate convolution identity, i.e., we 
assume that 4)*pt ^> as t for all e L^(?#). We will assign families of convolution kernels 
V'p and 4'p to pt such that 

/oo 
^/'p**pa(p)dp. (3.1) 

We assume that both families belong to L^(?^). Then we consider the following pair of transfor- 
mations. Let G L^(^^). Then we can assign to a two-parameter function Wcf) depending on 
points g ^'S and a scale parameter p G R+, the wavelet transform 

m ^ W^{p,g) ^(t>*M9) = / Hh)^^Jp{h-'g)d^l^{h) = (0,T;Vp), (3.2) 

and similarly we can invert this formula by convolution with ^Pp and integrating over the scale 
parameters 



W<P{p, •) * *pQ(p)dp = * / tjjp* *pa(p)dp. (3.3) 

Here and later on %o' as integral limit means that the integral has to be understood in an 
improper Riemann sense. 

This approach works for arbitrary approximate convolution identities pt. Of interest are in par- 
ticular those for which the operator 't^dtPt is positive. Then the corresponding Fourier coefhcients 
are positive matrices and the choice ^pp — ^p seems reasonable. We will later implement this 
general philosophy for the particular situation where pt is the heat kernel and where both families 
coincide. 
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Definition 3.1. Let ^ C # be co-finite. A family t ^ pt from C^{R+; L^{'^)) will be called 
diffusive approximate identity with respect to ^^_|_ if it satisfies 

(1) ||Pt(7r)|| < C uniform in tt G #and t G R+; 

(2) liiRt^o Pt{Tr) = I for all tt G 

(3) lim4^oc?5f(7r) = for all vr G 

(4) —dtptiTr) is a positive matrix for all t G M+ and vr G 

Remark 3.1. Conditions (1) and (2) of Definition 13. Il imply by Banach-Steinhaus that (/)*pt — 
as t — ?• in The norm || • || in (1) is the Euclidean matrix norm. Condition (1), (3) and 

(4) yield Pt\cf^ — ~ Jt^°° dpPpdp as L^(^^)-valued Riemann integral. 

Remark 3.2. If pt is a convolution semi-group, pt*Ps = Pt+s, the conditions are equivalent to 
the negativity of its generator L on (keriZ)^ — ^^^^^ "Htt- In particular the heat kernel pt as 

defined in (|2.11l) is an diffusive approximate identity with respect to = {vr G : > 0} = 
{tt G -H^ 7^ C}. 

3.2. Diffusive wavelets. For the following we fix a diffusive approximate identity and denote 

Llm - n^- (3.4) 

Wavelet transforms will live on that subspace of square integrable functions. Furthermore, we 
will require the admissibility condition (j3.1[) for the projection of pt onto this subspace. In 
the situation of the heat kernel the space L}^{^) contains all square integrable functions with 
vanishing mean. Furthermore, we denote by (/>|^ the orthogonal projection of G L^{^) onto 

Lq{^) and its obvious extension to distributions on J^. 

Definition 3.2. Let pt be a diffusive approximate identity and a{p) > a given weight function. 
A family ipp G Lq{^) is called diffusive wavelet family, if it satisfies the admissibility condition 

/•OO 

Pt\^^= Ipp * ippa{p)dp. (3.5) 



Equation (j3.5p can be solved explicitly. Applying the group Fourier transform to both sides 
yields 



Ptiir) = VpW^;(vr)a(p)dp, Vtt G (3.6) 

and by differentiating both sides we obtain 

-dtPt{n)=M^Wp{7r)a{p). (3.7) 

Hence, up to a unitary matrix ?Ap(7r) is the positive square root of the matrix — 9tpt(7r). 

Example (heat wavelet families). Considering the particular example of heat kernel pt and re- 
sulting heat wavelet families, we can go a step further and write (13. 7p as 

i^piiT) ^ -^X^e-P^l/^rj^ip) (3.8) 

for any (fixed) choice of a family rj^^ip) G \J{dTr). This implies 

iPp = — i= V d^X^e~P^'/^ trace(7r(g)??^(p)), (3.9) 
\/a(p) ^ 
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The freedom to choose rjj^ip) should not be overrated here; since we want to have Tpp concentrated 
around the identity element of the group the most sensible choice will be rj^^ip) = I. 

The weight function a{p) can be used to normalise the family V'p- If we want to have a normalised 
L^-norm, Parseval identity yields an expression for the corresponding weight a(p) 

a{p) = ^ d^A^e-"^^ = = -A^pp(l) ^ p-^-^ '^'-^ (3.10) 

7re# 

based on l|777rl||f5 = trace I = c?7r- Note, that Pp(l) is just the heat trace on . Another notable 
possibility is to use a[p) = 1 and rjT^{p) — rj^^ independent of p such that ipp = Pp/2 * ipo for 
^0 = E^c?-A,trace(^*(g)r;,) G P'(^). 

Theorem 3.3. Let be a dijfusive wavelet family. Then the associated wavelet transfrom (j3.2D 
W(t>{p, ■) — (f) * Ipp is unitary, 

W : L^(^) ^ L^{R+ X ^, a{p)dp (E)dpc^). (3.11) 
Proof. Let (pi and 02 be two functions with vanishing mean. Then 



{w^i,w<P2)= // WMp,g)'^Mp.9Wwi9Mp)dp 



cj)i{h)'ipp{h ^g)(t)2{h')ipp{h' ^g)dp<^{h)dp<^{h')dp<^{g)a{p)dp 

POO 



= lim / / (pi{h)(p2{h') I I ipp{h g)-ipp{h' ^g)dp<f{g)a{p)dpdp<^{h)dpc^{h') 



^}^^]] Mh)Mh') {Pt{h-^h') - 1) dp.^{h')dp.^{h) 
= {(Pi.h) 

based on cp * pt ^ (p in Lq{'^) as t — ?> and the statement follows. □ 

Remark 3.3. The functions ipp and Tgippi are not orthogonal in general. A simple calculation 
yields for heat wavelet families the identity 

(V'p,r,Vp') = , , \ , E rf^A2e-^'^x^(ff) = - , , A^p^(g). (3.12) 

4. Diffusive wavelets on homogeneous spaces /Jif 

We can use the wavelet transform on to define corresponding transforms on homogeneous 
spaces of ^. A particular example we have in mind are spheres §" as homogeneous spaces of 
S0(n + 1). 

Let for the following A" be a homogeneous space of i.e., we assume that there is a transitive 
group action '^^ x X ^ X denoted as {g^x) ^ g ■ x and satisfying gi ■ g2 ■ x = (5132) • x. If we fix 
a base-point xq X functions on X have a canonic interpretation as functions on ^ via 

4>ig)^H9-xo). (4.1) 

Note, that this identification depends on the chosen base-point in an essential way. Functions 
obtained in that way are constant on cosets gJff for -Jif = l^^o — {h ^ 'i^ : h ■ xq = xq} the 
stabiliser of xq in 5^. 
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We can define a {^-invariant measure dx on X by setting 

c^{x)dx^ I 0(g)dM<^(g) (4.2) 



and ask for a wavelet transform within L {X) 



Lemma 4.1. Equation (|4.2I) is independent of the choice of the base-point xq and defines a 
-invariant Radon measure on the homogeneous space X. 

There are two possible ways to define the wavelet transform. The first one is to apply the 
previously defined transform to the lifted function on the group the second one is to define a 
genuine transform associated directly to the homogeneous space. We will discuss both approaches 
in the next subsections. For simplicity we will first treat the simpler situation where wavelets on 
the group are expanded in terms of characters. 

4.1. The naive way. We now apply the wavelet transform to the lifted function for some 
4> & LP'{X). This defines a function on IR+ x via 

W^{p,g)^ f i{h)i;p{h-^g)dii^{h)^ f ^{h ■ xo)^p{h-' g)d^l^ih). (4.3) 



This definition has one obvious drawback; we would prefer to have a transform living on IR+ x X 
instead of R-^ x ^. 

Example (heat wavelet families). Heat wavelet families can be defined in terms of characters, i.e. 

with a (fixed) choice of a family rj^rip) G C, |?77r(/')| = 1- Then ipp is a central function on 
i.e. it is constant on conjugacy classes. Then the naively defined wavelet transform lives on the 
homogeneous space X, 

Wx^ip^x) - W4>{p,g) = / 4>{h- xo)iPp{h~^g)dtx^{h) 



= / (l){h-x)ijpih)dfi^{h) (4.5) 

where g is such that x — g ■ xq. 

The same reasoning works as long as the wavelet family on J# consists of central functions. We 
will refer to this transform as (naive) zonal wavelet transform on X. We can write (j4.5p as 
integral over X. Indeed, averaging formula (|4.5p over yields 



Wx4'ip,x)= I (j){h-XQ) / ^ppig ^hh')dpjp{h')dfii^{h) 



^ / tpih ■ xo)Vx^p{g ■ xo)dfi<:f{h) 



IX 

with 

Vx^p{h-xo)= I ^p{hh')dpMh') (4.7) 
J 

the ^-average of V'p (which is constant on cosets hM') seen as function on the homogeneous 
space X with fixed base point xq. The measure dp^ is the (normalised) Haar measure on M' . 
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Corollary 4.2. Let tpp be a diffusive wavelet family and assume that tpp are central functions on 
'S . Then the associated zonal wavelet transform Wx is unitary Lq(X) — > L^(R+ xX, a(p)dp'^dx) . 

Proof. The calculation is reduced to the one on Indeed, 



{Wx(bi,Wx(j)2) = / / WxMP,xWxMP,x)dxa{p)dp 
J^a J X 

= {W4>i,W4>2) = (01,02) = (01,02) 

follows from Theorem ES] and (g^])- □ 

4.2. Zonal wavelet transforms. We want to have a closer look at the structure of our formulas. 
For this we recall some notation. 

Definition 4.3. A function on X is zonal with respect to the base point xq if it is invariant 
under the action of the stabiliser of a;o. 

This is the case for Pa'V'p true in general for any projection of a central function x on 
¥xx{h-x)=¥xx{hg-xo)^ [ x{hgh')dp^{h') 

x{gh'h)dp,^{h') = Fxx{x). (4.8) 

To avoid a too cumbersome notation, we will simply and write instead of for Fourier coeffi- 
cients of functions on X (meaning that we identify X for fixed base-point xq with 'i^/J^). 

Proposition 4.4. There exists a family of orthogonal projections J^{tt) G C'^vrXd^ depending on 
the stabiliser such that 

f^^{Tr) = ^(7r)0(7r). (4.9) 

Proof. A short calculation yields 



J^(7r) = / 7:{h)dpMh), (4.10) 

such that ,yf{TT),jr{TT) = ,yf{TT) aud ,yf*{n) = .^{-k). □ 

Corollary 4.5. A function G L^{X) is zonal if and only if J^{tt)(P{tt) — 0(7r)J^(7r) for all 
TT e This is true, if and only if there exists a function £ L^{X) with = 0. 

Since ^(tt) is an orthogonal projector in C^^ we can choose a basis of C*" in such a way 
that J^{tt) — b-diag(Ir^ , 0) is block-diagonal, — rankJ^(7r). Then, for S L^{X) Fourier 
coefficients satisfy (I — J^(7r))0(7r) = 0, which means, it is a matrix having non-zero entries in 
the first r^r rows only. Similarly, e L^{X) is zonal if its Fourier coefficients have entries only 
in the upper-left x block. 

Convolutions of functions defined on X should be related to multiplications of Fourier coeffi- 
cients. The underlying block structure of Fourier coefficients makes it advantageous to combine 
multiplications with taking adjoints. This leads to the introduction of several convolution like 
products for functions; we collect them in the following definition. 

Definition 4.6. Let 0,0 G L^{X). Then we define 
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(1) the group convolution 

^ * i^ix) = I ^{g- a;o)V(g-i • x)d^l^ig) e L^X); (4.11) 

(2) the •-product 



4> • ^{g) = / ^{x)i^{9-^ ■ x)Ax = ((/., t;^) e L\'^)- (4.12) 



(3) the zonal product 



cj>i^(x) = / c^{g-xMg-x)Aii^(g)(,L\X). (4.13) 



In order to justify the above definitions and to relate them to Fourier coefficients we collect some 
properties of these convolution like products in the following proposition. 

Proposition 4.7. Let (l),ip e L^{X). 

(1) (p * ip = (f) * and thus (j) * tp{Tr) = '0(7r)(^(7r). 

(2) (j) u tp = (f) * and thus p • ip{n) = 'iP*{'k)P{tt). 

(3) // Tp is zonal with respect to xq then (p • ip is constant on cosets gJ^ and thus defines a 
function on X. 

(4) (piip is zonal with respect to Xq ■ 

(5) (piip — 4> * ip and 4>itp{T:) — iJj(7t)(P*{tt). 

Proof. (1) evident. (2) It suffices to consider the first formula which follows directly from (|4.2p . 
(3) We calculate 



p»Tp{gh)— / (p{x)ip{h ^g ^ • a;)da; = / (p{x)ip{g ^ ■ x)dx — p • 'p{g) 
J X J X 

based on the zonality of ip. (4) For h £ = 'iS^o we obtain 



<pi-p){h-x)= \ p{g ■ xo)ip{gh ■ x)dfj.^{g) = / 4>{gh ^ ■ xo)ip{g ■ x)dfj.':^{g) = piipix) 

J'S J'S 

by the right invariance of /i^. (5) It suffices to show the first formula. Substituting g~^ for g 
and using (j4.2l) implies the assertion. □ 



Corollary 4.8. Let (pj'ip ^ L^{X). If ip is zonal then <p * tp = p • ip. Similarly, if p is zonal then 

(p ^ 'p! — (p i xp . 

Corollary 4.9. The (naive) zonal wavelet transform is given by a •-product 

Wx4>{p,x)=P^¥xi>p (4.14) 

with inversion 

0(x)=/ WxPip,-)»^x^Jpaip)dp. (4.15) 



The •-product of zonal functions is not associative. However, a short computation (e.g., on the 
Fourier side) yields [p • ip) • x = • (x * V')- Therefore, the admissibility condition implies or 
rather reads as 

/oo 
¥xi^p*¥xiPpa{p)dp. (4.16) 

In fact, this can be used as starting point to define zonal wavelet transforms on X . 
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Definition 4.10. Let pt be a diffusive approximate identity and let a{p) > be a given weight 
function. A family ijjp € L'^i'^) is called zonal diffusive wavelet family if 

(1) ipp is zonal with respect to xq, 

(2) the admissibility condition 

/oo 
i^p * i;p{x) a{p)dp (4.17) 

is satisfied, where 

Pt^ig-xo)^ / Ptihigh2)dfi^{hi)dp,^{h2) (4.18) 
is the zonal average of pt on X. 
We associate to this family the wavelet transform Wx4>{Pi x) = (p* ippix) with inversion given as 

Wxc^{pr)»i^pa{p)dp (4.19) 



for all functions £ L/^{X). 

Theorem 4.11. The zonal diffusive wavelet transform is unitary Lq{X) — > L^(R+ x X, a{p)dp<S) 
dx). 

Proof. This can be shown along similar lines to the proof of Theorem l3.3l Alternatively, one can 
use Plancherel's theorem and show that the corresponding transform in Fourier space is unitary. 
We follow that line of thought and write 

(</'i,</'2)= ^ trace((^2(^)'^i(^)) = X/ ^^tt trace((^2(^)-^(^)'^i(^)) 
^ d„trace(02(7^)V'p(7r)V'*(7r)0i(7r))a(/?)dp 

= (#A"^l,^^</'2) 

based on (liTTl) in the form of 4>p{TT)^/j*{TT)a{p)dp = ,3^{tt) for tt □ 

4.3. Nonzonal wavelets. Now we will consider the situation where ■(/'p G L'^i'^) is an arbitrary 
(non-zonal) family. Then 

W<P{p, g)^(b» Vp(5) = Vp) (4.20) 

lives on W rather than on X. We aim for an inversion formula of the kind 

/oo 
Wcf>{p,-)*^p{9)a{p)dp (4.21) 

with a second family G L^{X). By a short computation we see that {(j) • * x — 4> ' (x* "0) 

for (j),ip,x G i^('^) and this motivates us to define 

Definition 4.12. Let pt be a diffusive approximate identity and a{p) > be a given weight 
function. A family ipp £ L'^{X) is called (non-zonal) diffusive wavelet family if the admissibility 
condition 

/oo 
iPpi^p{x)a{p)dp (4.22) 

is satisfied. 
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We associate to this family the wavelet transform Wx4>{p,g) = 4> • ippig) with inverse given as 

/oo 
l^x<l>{p,-)*i'pa{p)dp (4.23) 

for all (j) e Ll{X). 

Remark 4.1. If is zonal, i'p*ipp = ^p* ipp and Wx<j){p, ■)*V'p = '^x4'{p^ O'V'p by Corollarv l4.8l 
Hence Definition 14. 121 extends Definition 14.101 

The proof given to Theorem 14.111 transfers to the non-zonal situation. Hence, we obtain 

Corollary 4.13. The (non-zonal) diffusive wavelet transform is unitary Lq(A') — S> i^(R+ x 
a{p)dp ® d^gp). 

Remark 4.2. We want to emphasise the compatibility of the wavelet transforms on X and on . 
Therefore, we look at the transformations in Fourier space. The admissibility equation reads in 
both cases as an integral of iPp{tt)iP*{tt). If ipp is a wavelet on X the result are Fourier coefficients 
of a zonal function. 

The wavelet transform on as well as on X can be written in Fourier space as 

W^ip,7T)=rp{7r)4>{7r). (4.24) 

On X this gives in general a Fourier coefficient for a function on (i.e., a full matrix). Wcf) is the 
Fourier coefficient of a function on X if and only if -tpp is Fourier coefficient of a zonal function. 
In both cases the reconstruction is given as integral 

/oo poo 
V'p(7r)^(p,7r)a(p)dp= / i^p{7T)rp{n)a{p)dpf{7r). (4.25) 

Example (heat wavelet families). We can describe all admissible non-zonal heat wavelet families 
by looking at their Fourier coefficients. Formula (|4.22p implies V'p(^) — for Att = and 

Xljf{7T)e~^"p - V-p(^)7A;(^)a(p), A, > 0. (4.26) 

We assume that we have chosen the basis in C^" in such a way that J^(7r) ~ b-diag(Ir^ , 0), 
r-rr = rank J^(7r). Then, ipp{7r) has non-zero entries in the first r^r rows only and condition (j4.26p 
implies for these entries that (i) the rows are pair-wise orthogonal 



^{M^)UiM^))jk ^ 0, i^j, i,j^l,...,r^ (4.27) 

k=l 

and (ii) of prescribed length 

E = ^^A^e-^^^", * = 1, . . . , r.. (4.28) 

Of particular interest are situations where G {0, 1}. Then (I4.27P disappears and (|4.28p reduces 
to one condition. This corresponds to the assumptions used in jS] and [3]. 

We want to remind the reader of the following definition, which is particularly helpful when 
considering concrete geometric situation as we do in Sections 15.21 and 15.31 

Definition 4.14. The subgroup Jif is called massive, if lankJ^iir) < 1 for all tt e Further- 
more, an irreducible representation tt G is called class-1 representation with respect to the 
subgroup J^, if rank .^(tt) > 1. 
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5. Examples 



5.1. Tori and functions. First, we consider the simplest possible situation. Let 'S = T = 
U(l) C C be the set of unimodular complex numbers. Then T = Z and to each k (ilj corresponds 
an irreducible representation z ^ . The corresponding Fourier series are just Laurent series 
or the usual Fourier series if we write z — exp(27ri0). The Laplacian on T corresponds in this 
notation to dg and has eigenvalues — 47r^fc^. 

In this notation, the heat kernel on T is given by 



p^{z) = e"*"''''*^'' = 1 + 2 ^ g-4^=fc=t cos(2fc7r0) = M^9, e-^-'*) (5.1) 



kei fc=i 
in terms of Jacobi's z?3-function, cf. 26, Chapter XXI]. Similarly, one obtains 



y^i/,p(z) = -27ri^fce-2^''='''z'= = 47r ^ fce'^'^''''^ sin(2fc7r6l) = -deMT^^^e"^""' ") (5-2) 



(choosing rjk — — isignfc to simplify the notation). The corresponding wavelet transform of a 
function / G L^[0, 1] — L^(T) with normalisation a{p) = 1 is 



pi 

= / /'(0-T)i^3(^T,e-2"''')dr 
Jo 



(5.3) 



with inversion formula 



</>(r)dT- / / W4>{p,9 ~T)drM-^T-,e-^'''P)dTdp. (5.4) 

The wavelet transform W(f>{p, 9) describes for small p the 'high-frequency part' of ip localised 
near the point 9. 

This t9-transformation generalises to higher dimensional tori in the obvious way. We do not give 
further formute, but conclude this example with some pictures of the family ■0p for different p 
depicted in Figure [TJ 



5.2. The group S'^. The sphere can be viewed as group of unit quaternions, S'^ C H, or 
alternatively, as group of unitary 2x2 matrices SU (2) (see, e.g., [53] for that approach). We 
decide for quaternions x = xq + ixi + ix2 + kx^ G H with Xi G M, i^ = = = ijk = —1. The 
unitary dual S'^ can be identified with Nq ; for each n G No there is a representation 

tn : §^ ^ U(n„), tnix)pn = Pnix''^-), (5.5) 

on the space n„ C C[xo, - ■ ■ ,2^3] of all polynomials in four variables which are homogenenous 
of degree n and L^(S'^) perpendicular to n„', n' < n. As IJn n„ is dense in L^(§^) we obtain 
all representations in this way. Furthermore, n„ splits into dn minimal invariant subspaces of 
dimension dn corresponding to the splitting of i„ into irreducible components. 

An orthonormal basis of n„ is given by spherical harmonics, which occur as matrix coefhcients 
of our representation {7^*-' (a;) : i,j — 1, . . . dn}, dn = n + 1. As the maximal torus of S'^ is of 
dimension 1 (and we may fix it as n C = {x G S'^ : 2:2 = X3 = 0}) class functions are functions 
of just one variable a;o = Sex. This is seen easily from y~^xy = Sc x + {Yec x)y . Furthermore, 
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Figure 1. The toroidal family 1^3(776', e for -| < < | and scale param- 

eters p e {0.005,0.01,0.015,0.025,0.04,0.1}. 



such functions are essentially unique within n„; they are multiples of the Gegenbauer polynomial 
C^(xo). Hence, using C}^{1) = n + I — dn we obtain 

d„ 

trace i„(a;) = ^^X^^ {x) = C^^ixo), xq = Sex. (5.6) 

We normalise the measure on S"^ to be the Haar measure (which changes the Laplacian from the 
usual one ). Then Xl = (27r2)2/3„(„ + 2) and we obtain the corresponding heat kernel on S as 

00 

PtW = Xl('^ + l)^~^"*^"(^o)' a;o = Sca;, (5.7) 

n=0 

such that the corresponding heat wavelet family is given as 

00 

^p{x) = ^== V(n + l)A„e-^"''/2ci(xo), xq = Sex. (5.8) 

V "(P) n=l 

Figure [2] depicts ijp{x) on S'^ n C for different p. 

Some spherical 3-manifolds can be treated similarly to S'^ in the previous section. Let for this 
F be a finite subgroup of S'^ C H and X = S'^/F. Interpreting all integrals over F as averages, 
we can retrace the considerations of Sections 14. II to A zonal heat wavelet family on S^/T is 
therefore given by 



M^) = ^^Y.{n + l)Ke-^'^^/' I 1 ^ci(Scx7) I 



(5.9) 



Main ingredient to understand zonal and non-zonal transforms on the Fourier level is the family 
of matrices 

intn)h^^,J2T2'{i). (5.10) 



|r| 
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Figure 2. The constructed heat wavelet family on S'^. For the figure we have 
chosen a{p) - p-^ and p € {O.OOf , 0.003, 0.005, 0.007, O.Of , 0.02, 0.03}. 

As example we consider the lens spaces L{p,l). Then r=:{ti;eS'^nC:iuj'° = f}~ Z/pZ and a 
simple calculation (based on [SJ 10.8(17)]) yields 

1 ^/ \ ^/ \ 1 ^1 / 2fc7rs^ 1 2(n — 2m)fc7r 
rankr(t„) — trace r(t„) — - > C„(cos j = - > > cos 

^ k=l ^ ^ k=l m=0 ^ 

= #{m e {0, . . . , n} : n - 2m EE 0{p)}. (5.11) 

In particular for the projective space RP^ = L(2, 1) it follows that r(t„) = I for n even and 
r(in) = for n odd. The zonal heat wavelet transform on projective space is thus obtained by 
ignoring all odd Fourier coefficients in ([57 



5.3. SO(?T. + 1) and Part of the computations for this example are taken from We will 
just recall some notation, for more details see or [24|. As usual we denote by yi G C°°(§"), 
fceNo,j = l,..., dk{n) the orthonormal system of spherical harmonics on S" (normalised with 
respect to Lebesgue measure on §"). Their span is invariant under the action of SO(n + 1) and 
therefore it makes sense to define the Wigner polynomials 7^"'(A), A € SO(n + 1), 

dfc(n) . 

yi{A-'0=Y.V.'{A)yi{i), Ti'{A)= yi{A-^oyl{m (5.12) 

as matrix coefficients of the corresponding representation. The matrix-valued function A H> 
Tk{A) forms an irreducible unitary class-1 representation of SO(n + 1) on §" ~ SO(n + 1)/S0(ri). 

Remark 5.1. The Wigner polynomials do not give all irreducible unitary representations of 
SO(n + 1). As we are mainly interested in §" here, we do not need to care about that. Lifts of 
functions on §" to SO(n + 1) have vanishing Fourier coefficients on all other representations. 

Lemma 5.1 ([24, Chapter IX. 2]). S0(7i) is a massive subgroup o/SO(n + l). Furthermore, the 
family 71-, k G Nq, gives up to equivalence all class-1 representations o/ SO(n + 1) with respect 
to SO(n). 
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For the following we fix the 'north pole' of S". Then it follows that the set of zonal 
spherical harmonics (w.r.t. ^o) is one-dimensional and spanned by the Gegenbauer polynomial 
^(n i)/2^^T^^^ This implies several interesting formulae, in particular 

Lemma 5.2 (Addition theorem). For all ^,77 G S" and fc e No 

C<'— >"(«T,) _ IS.., 



Sketch of proof. It suffices to check that the right hand side is zonal with respect to ^, which 
follows from (j5.12p with A in the stabiliser of ^ and exchanging orders of summation. Then, in 
order to find the constants it suffices to choose £, = rj and integrate both sides over □ 

Lemma 5.3 (Zonal averaging). For all G A; G Nq and i = 1, . . . , dk{n) 



("-l)/2/fT 



Sketch of proof. As the left hand side is after averaging zonal with respect to the north pole £^0, 
it is again just a matter of determining constants. For this we set C = Co- 

Lemma 5.4 (Funk-Hecke formula). Let f : [—1,1] — > C 6e continuous. Then for all i = 
l,...,dk{n) 

/(g^^)K(Odg^K(^) J_i)J . / f{t)ci''-'^/'{t){l-t^r/^-'dt. (5.15) 

Sketch of proof . As /(^^f?) is zonal with respect to 77, we can average the integrand by replacing 
77 with Br] and integrating over B from the stabiliser of 77. Using the invariance of the integral 
on S" this can be transferred to an average over yi and therefore (|5.14p reduces the integral 
to an integral over fi^^rfjC^ ^^^^(C^*?)- It remains to substitute t — ^^77 and the statement 
follows. □ 

5.3.1. Zonal heat wavelets on S". At first we want to give a formula for zonal wavelets. Again we 
fix ^0 to be the north pole of S" and set ^ = Afo- Then the S0(7i)-averaged character trace 7fc(A) 
has to be zonal w.r.to and is therefore given by 

(P..tracer.)(0 = ^^^^ trace (I) = ^^^ct'"\e 

= ^-^^Cr^''\ei.) (5.16) 

based on d,{n) ^ - ("+^2^ and ct'-'^'\l) = ("+^^). 

Hence, a family of real zonal heat wavelets on §" is given by 

= ^ E ^^i^A.e-^-/^Cr-)/^C^Co), (5.17) 

where — is the corresponding eigenvalue of the Laplacian on SO(77 + 1). 
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Remark 5.2. In fact, any other sequence of positive numbers Xk oo will give a suitable wavelet 
transform with similar properties. In this situation we would replace the heat semigroup with an 
arbitrary semigroup generated by a non-negative bi-invariant operator (given by the sequence of 
eigenvalues A^, on its eigenspaces span{7^-' : 1 < i,j < dk{n)}). 

5.3.2. Nonzonal wavelets on To shorten notation let — SO(n). We first calculate the 
projection Jif{k), i.e., 

7sO(n) \JS" JsO{n) ) .. 



X^k "^"1 J 



l§"l 

dfc(n) 



(K(6)3^^(eo))^. (5.18) 



based on (|5.14p . Funk-Hecke formula and the normalisation of Gegenbauer polynomials, cf. [SJ 
11.4(17)], 

^1 r^(n-l)/2/, Ni2 

|g„-i| / p(n-i)/2^^^j2^^ _ t^yl^-^At = 1^L^^|§«|. (5.19) 

In order to simplify notation, we assume that the basis of spherical harmonics 3^^(C) is chosen 
in such a way that yl{£.o) = V'^fe(n)/|§"| and 3^^(Co) for aU i > 1. Then Jf'(fc) = 
diag(l, 0, . . . , 0). In particular this means 



dk{n) 



ym^{p^-v'm= f T,:\AB)dB = T^\A), ^^a^o- (5.20) 

JSO(n) 



In order to specify a non-zonal wavelet family we have to specify the corresponding Fourier 
coefficients tpp{k) for each class-1 representation Tfe. Due to the considerations in Section l4!3l 
we know that ippik) has entries only in the first row and these entries satisfy (|4.28p . To single 
out an interesting family of non-zonal wavelets we follow [3] and introduce an weight vector 
w{k) € cdfc(ri)^ |K;(fc)p = J2i = 1; define the Fourier coefficients in terms of w(fc) as 

4>p{k) = —l=XkC-^"-p/^ek(E)w{k), (5.21) 

where Ck = (1,0,..., 0)^ e cMn)^ _xl is the sequence of eigenvalues of the Laplacian on 
SO(n -I- 1) and a{p) a suitable weight function. The choice w{k) = Ck corresponds to the zonal 
wavelets constructed in Section [5.3.11 

For n = 2, dk(2) = 2k + l, we get therefore for any choice of weight vectors w{k) £ C^'^"'"^ of unit 
length a family of non-zonal wavelets, 

oo k 

MO = ^= E(2* + E M>')ykiO, (5.22) 

Va(p) fc=o i=-k 
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Figure 3. An example of a non-zonal heat wavelet on Scale parameter 
is p = 0.15, weight-vector wik) = {2k + . . . , 1)^ e C^'^+i is cqui- 

distributed and weight function chosen as a{p) = p^^ . 



where the spherical harmonics on § are denoted in the more conventional way as 



yliO = \r-^Cl'\tl (5.23) 



22-bV2r(fc+|j-| + i) ^^-b-l M e , .7- k,...,k, (5.24) 

in the special spherical coordinates t = S^q ^ describing the sine of latitude and 9 G [—t^, tt] any 
convention of longitude on S'^. Figure |3] depicts an example of a non-zonal wavelet on S'^. 

The associated non-zonal wavelet transform of a function / e L^(S^) lives on K+ x S0(3). We fix 
the north pole on §^ as point with co-ordinates (0, 0, 1)^ and introduce Euler angles a G [— tt, tt], 
/3 G [—n/2, 7r/2], 7 G [— tt, tt] as coordinates on S0(3), such that 

(cos 7 sin 7 0\/l ^\ / '^'^^ ^ 

— sin 7 cos 7 cos /3 sin /3 — sin a cos a 
1/ yO -sin/3 cos/3/ \ 1^ 

= R^S/bRc. (5.25) 
Then we can write the wavelet transform as 



#'/(p,a,/3,7) = /-^p= / /(0V'p(^-i(a,/3,7)0de 

JS2 



/ f{R^SpOMR-am- (5-26) 

JS2 
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Hence, /3,7 specify the point under consideration (the preimage of the north pole (0, 0, 1)^ under 
the rotation R^Sp), while a specifies directions (related to rotations of the wavelet ipp around 
the fixed north pole (0, 0, 1)^). 

6. Concluding remarks 

1. We defined wavelet transforms for square integrable functions on Lie groups 5^ and their 
homogeneous spaces X. Since heat wavelet families are smooth, all considerations generalise 
immediately to distributions u G 2?'(^^) or it e 'D'{X) and give rise to a smooth function Wu{p, g) 
decaying exponentially as p — > oo. 

Furthermore, the approximation properties of the heat kernel are directly linked to the inversion 
formula. Hence these are valid in a variety of function spaces, e.g., for any (f) 6 L'^i'^)^ P £ [l^ 

/•oo 

= lim/ #'0(p,-) * Vp(-)a(/o)dp inLP(^). (6.1) 

t-s-O 

For arbitrary diffusive approximate identities in the sense of Definition IS . 1 l onlv the L^-convergence 
follows. The transform is still well-defined on 2?'(^^) or 'D'{X). 

2. Regularity properties of functions can be expressed in terms of their wavelet transform. To 
be more concrete on this, we recall that (p G L^(^#) belongs to the Sobolev space iJ'*(^^) if the 
Fourier coefficients satisfy 

mh = J2dAKruM\\%s<^- (6.2) 

We say (p belongs locally in g G to H'^, if there exists a cut-off function x £ C°°(^^) with 
g G suppx and x(f> G H^{^). The latter property can be related to the asymptotic properties 
of the wavelet transform W(f>{p, g) evaluated in g as p — )> 0. Non-zonal wavelet transforms on 
homogeneous spaces allow the same consideration micro-localised to directions. 

3. It seems interesting to ask for discretisations of these wavelet transforms. As W4){p,g) is 
smooth in p > and 5, we can evaluate it in discrete points. If we choose them in an appropriate 
way the inversion formula imply that the corresponding wavelets form a frame in L^{^). We 
refer to a forthcoming publication concerning details of this. 

4. Spherical wavelet transforms are of interest for applications in geophysics, crystallography 
and medical imaging; see, e.g., [13] or [TH] and references cited therein. A particularly inter- 
esting situation arises for homogeneous spaces S0(3)/r for diskrete rotation groups F < S0(3) 
representing crystal symmetries. 
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